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We have studied a hybrid system of a one-dimensional coupled-cavity waveguide with a two- 
level system inside, which subject to a external periodical field. Using the extended Hilbert space 
formalism, the time-dependent Hamiltonian is reduced into an equivalent time-independent one. Via 
computing the Floquet-Green's function, the Zeno-anti-Zeno crossover is controlled by the driven 
intensity and frequency, and the detuning between the cavity and the two-level system. 
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I. INTRODUCTION 



Inspired by modern microfabrication technology in 



photonic crystalsjl|; 
ducting devices 



|, optical microcavitiesp, supercon- 
and the realization of the quantum 
regime in the interaction of atomic-like structures and 
quantized electromagnetic modes in those system [IoHl4) . 
considerable attention has been attracted to a family of 
models for coupled arrays of atom-cavity systems, since 
they offer a fascinating combination of condensed matter 
physics and quantum optics. Typically, these models are 
formed by an array of cavities with each cavity containing 
one or more atoms [1 51] . where photons hop between the 
cavities. Recently, due to the potential use for building 
a quantum switch for routing single-photons in quantum 
network, systems with one or two atoms inside the ar- 
ray of cavities have been extensively studied to revealed 
the intriguing features of photon transport in low dimen- 
sional environments fla - 1201 ] . It is found that the switch 
is formed by the interference between the spontaneous 
emission from atoms and the propagating modes in the 
one-dimensional (ID) continuum. Therefore, the search 
for a controllable switch is to finding a way to change the 
spontaneous emission of atoms. 

The spontaneous emission results from the inevitable 
interaction of the atomic system with external influences. 
Decoherence of a quantum system, with a variety of cou- 
plings to a reservoir, has been investigated extensively in 
theory j2l| - [23| . Basically, there are three ways to suppress 
or modify the rate of quantum transitions in a system. 
One is to engineer the state of the reservoir, as well as 
the form of the system-reservoir coupling [24H26j . Ob- 
viously, the coupled-cavity waveguide (CRW) meets this 
condition for its advantage of addressability of individual 
sites, extremely high controllability, and the great degree 
of flexibility in their geometric design. Two is to involve 
the quantum inference between multiple transition path- 
ways of internal states, for example, the electromagnet- 
ically induced transparency technique p7l - l29| . The third 
way is to applies a succession of short and strong pulses, 
or measurement to the quantum system [30U35| . 

Atoms have been particularly interested for acting as a 
quantum node in the extended communication networks 



and scalable computational devices, specially artificial 
atoms. Atoms in a time-varying field has been inves- 
tigated long time before. However, the multiphoton res- 
onance and quantum interference have been experimen- 
tally demonstrated in a strongly driven artificial atom [36| 
until recently, which is important to superconducting ap- 
proach of quantum computation, for instance, decreasing 
the time required for each gate operation |37H40| . 

Floquet theory is developed by Floquet centuries 
ago, it is a theory about the solutions of linear differ- 
ential equations with periodic coefficients [4lj|. Later, 
it is applied to the two-level system (TLS) with the 
time-dependent problem, which is discussed by Aut- 
ler and Townes(42l|. And then such kind of the pe- 
riodic time-dependent problem is reformulated as an 
equivalent time-independent infinite-dimensional Flo- 
quet matrix [43j. which is done by introducing the com- 
posite Hilbert space of square integrable and time- 
periodic wave functions. Now Floquet formalism is 
used as a theoretical tool to investigate time depen- 
dent phenomena [44} . In this paper, we study a two- 
level system (e.g. a flux qubit) interacting with a cav- 
ity which together with other cavities constructs a one- 
dimensional (ID) coupled-cavity waveguide (CRW). The 
CRW is modeled as a linear chain of sites with the 
nearest-neighbor interaction. Obviously, the dynamic of 
the TLS is irreversible, i.e. once the TLS is initially in 
its excited state, it never returns to its initial state spon- 
taneously. Therefore, the CRW with a TLS inside is a 
typical system with a discrete state coupled to continua 
of states, which means the TLS is subject to decay. In 
order to control the decay rate of the TLS, an external pe- 
riodical field is applied to the TLS through diagonal cou- 
pling. It is well known that periodic coherent pulses can 
either inhibit or accelerate the decay into its reservoir, 
however, such modification is based on off-diagonal time- 
dependent couplings between the TLS and the external 
field. In a atomic experiment, the off-diagonal couplings 
is caused by a transverse field perpendicular to the po- 
larizing field, which is used to polarize the atomic spins. 
Then the spins of atoms and photon are aligned. How- 
ever, diagonal couplings between the driving field and the 
TLS means that the driving field is parallel to the polar- 
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FIG. 1. (Color online) Schematic illustration of the model, 
where a two level system is inside a ID coupled cavity waveg- 
uide. The two level system is driven by external forces that 
are periodic in time with period v. 

izing field (45j. Here, the spontaneous emission of a TLS 
is studied, the Zeno-anti-Zeno crossover is investigated 
via Floquet formulation. 

The paper is organized as follows. To set our consider- 
ation, in Sec. UH we state the Hamiltonian for ID CRW 
with a TLS inside one of the cavity, in Sec. IIII1 we first 
give a brief outline of the Floquet theorem, then present 
the Floquet representation of the system we consider. 
In Sec. ??, we calculate the amplitude for the TLS in 
its excited state in the one quantum subspace, and the 
quantum Zeno and anti-Zeno crossover is investigated. 
Conclusions are summarized at the end of the paper. 

II. DESCRIPTION OF THE MODEL. 

Consider a TLS subject to an harmonically transverse 
driving with frequency v and intensity A. The lower and 
upper eigenstates of the two-level system are described 
by notation |e) and \g) respectively, which is separated 
by energy f2 without external fields. Such kind of TLSs 
can been experimentally realized using superconducting 
circuits [3(1 [45[. The Hamiltonian describing the TLS is 

H A = {tt + Acosvt)a z /2, (1) 

where a z describes the atomic inversion, + A cos vt is 
the energy splitting. The TLS interacts with a quantized 
electromagnetic field of a ID waveguide, which is con- 
structed by coupling of cavities in an array. Due to the 
overlap of the spatial profile of the cavity modes, pho- 
ton hops between neighbouring cavity. Introducing the 
creation and annihilation operators of the cavity modes, 
and aj , the CRW is modeled as a linear chain of sites 
with the nearest-neighbor interaction. The Hamiltonian 
of the ID CRW can be written as 

Hw = £ ^c a \ a i ~ £ £ ( a j+i a J + h - c ) ' ( 2 ) 

j 3 

where £ is the hopping energy between adjacent cavities 
and 10 c is the eigenfrequency of each cavity. The TLS is 



inside one of the cavity in the ID CRW, which is labeled 
as the zeroth cavity. In the rotating wave approximation 
(RWA), the interaction between the TLS and the zeroth 
cavity is described by a Jaynes-Cummings Hamiltonian 

Hi = g (a-al + /i.e.) . (3) 

The operators a+ and ct_ are the usual raising and low- 
ering operators of the TLS. The Hamiltonian governing 
the system is 

H = H A + H w + Hi. (4) 
By employing the Fourier transformation 

the Hamiltonian H in the Hilbert space of configuration 
is now given by in the momentum space 

h = y^ £fc4°* + \ i n + A cos (^)i a * ( 6 ) 

k 

where the dispersion relation 

Sk = oj c — 2£ cos k (7) 

describes an energy band of width 4£ (the lattice con- 
stant is assumed to be unity). The second-quantization 
operators afc/a[- annihilate/create one photon in the kth 
mode of the ID CRW. The periodical boundary condition 
is used to obtain the Hamiltonian H in Eq.©. 

It can be found that the number of quanta, which is de- 
fined by operator j\f = J^ fe aj^a^ + a z , is conserved in this 
system, i.e. j\f commutes with Hamiltonian H . There- 
fore, if we have one quantum in the initial state, the state 
vector evolves restrictedly in the one-quantum space. We 
introduce the states |fc) = a\ \0g) which describes that 
there are one excitation in the fcth mode of the CRW 
while the TLS stay in its ground state, and |e) — |0e) 
which denotes that the TLS has been flipped to its ex- 
cited state while the CRW is in the vacuum state. The 
orthonormal basis set {\k) , |e)} spans the one-quantum 
space. Therefore Hamiltonian H is rewritted as 

ff=i[fi + AcoB(i/f)] U){e\-J2\k)m^j (8) 

+£ e *l*X*l + -^£(l s > < e "i + /i - c -) 

fc v fc 

in the one-quantum subspace. Hamiltonian H in Eq.® 
describes a single discrete state is coupled to a continuum 
when the periodical-driven field is absent, which means 
the discrete state is subject to decay. Consequently, the 
excited state of the TLS is an unstable state. 



III. FLOQUET FORMULATION FOR DR 
ATOM INSIDE A ID WAVEGUIDE. 



Since the external driving field is strictly per 
time, Hamiltonian in Eq.® is a periodic fun< 
time, i.e. H (t) = H (t + T) with T = 2n/v be 
period. To study our problem, it is necessary 
sider solutions of Schrodinger equation with a t 
riodic Hamiltonian. In this section, the theoretic 
that we used is the Floquet representation of qr 
mechanical systems. A brief outline of this met] 
been given in appendix. 

We now employ the Floquet-state nomencla 
| an) = | a) (8 |ra), where n is the Fourier index 
from — oo to oo, a = e, k is the system index, 
deal with Hamiltonian Ha in Eq.(fT|). In the 
space, the time parameter is another degree of 1 
of the system. Hence, the periodical function co 
an operator with the expression 

cos (vt) = i ^2 (\ n + !) ( n l + h - c -) ■ 

n 

In the one quantum subspace, Hamiltonian Ha — 
be separated into two segments H e + H g with 
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It can be diagonalized by the following transform 

| era) = ^2 J n -m (-x/ 2 ) l#m) (Ha) 

771- 

|fcn> = ^ J„_ ro (x/2) |fc0 m ) (lib) 

?n 

where J„ (x) is the Bessel function of the first kind and 
X = A/v. In terms of states in Eq.tfTTj) we reduce the 
solution of a periodic time-dependent Hamiltonian H to 
the problem of diagonalizing the time-independent Flo- 
quet Hamiltonian Hp = Hq + H\ with 



kin 



Hl= J2 gJm 7g (x) (l^ m )(e^l+fe.c.). (12b) 
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Here, £^ (i?^) is the eigenvalue of H e (H g ) with 

£4 = fi/2 + mi/, (13a) 
Efr =s k - + mi/. (13b) 
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FIG. 2. (Color online) Schematic Quasienergy diagram of the 
Hamiltonian H e (a) and H g (b) in Eq.fTO]) under the condition 
2£ < v. 



The eigenvectors of Hamiltonian H e and H g are coupled 
via the nonzero coupling strength g. Figure ?? shows 
the energy diagram of Hamiltonian H e and H g in Eq. fTU)) 
when 2£ < v. When the eigenvalues and are close 
to each other, i.e. E^ « E^, Floquet Hamiltonian Hp is 
reduced to the following form 



H R = H 



BO 



H 



(14) 



where 



H R0 = E e m \e4> m ) 



(e4> m \ +J2 E n \ k ^n) (k<p n \ (15a) 



gJn-m (X) 



N 



\k(p n ) {e(f> m \ +h.c.) (15b) 



This approximation seems equivalent to the rotating 
wave approximation (RWA) which is traditionally used 
in the field of atomic field to neglect the counter-rotating 
term of the harmonic driving, however this approxima- 
tion is different from the RWA, which is valid only for 
amplitudes of the driving field small compared to the 
energy difference between the atomic states and breaks 
down in the strong field. 

The above discussion shows that if the TLS is initially 
in its excited state, it will emit a photon as a result of the 
interaction with the radiation field of the CRW, photons 
will gain or loss energy quantum fno due to the periodical 
modulation. Hence the state of the emitted photon is 
characterized by the set of energy E q = f2 — qhuj with 
q = n — m. Once E q equals to the energy of the 
CRW, photons go to the CRW. However, when ratio of 
the intensity to the frequency of the modulation is equal 
to the roots of the Bessel function J q (x) , this process is 
prevented due to the decoupling of the TLS and CRW as 
one can see in Eq. (fT5b ). When the modulation is absent, 
the index q vanishes (i.e. q — 0). The behavior of the 
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emitted photon is determined by whether ft is equal to 
£fe or not. Here, we only give a intuitively discussion, 
more details will present in the next section. 



IV. ZENO-ANTI-ZENO CROSSOVER 

Time evolution of standard quantum theory assumes 
two principles: the continuous unitary evolution without 
measurement, and the projective measurement. Quan- 
tum Zeno effect (QZE) or anti-Zeno effect (AZE) is a phe- 
nomenon related to projective measurements, which says 
that repeated observations prolong or shorten of a life- 
time of an unstable state. However, slowdown or speedup 
of the decay of an unstable state is also possible without 
measurements or observations [tH l48l| . 

We now investigate the lifetime of the periodically 
driven TLS in interaction with the CRW when the hop- 
ping energy and the driven frequency are chosen such 
that 2£ < v. An excited state of the TLS in one quantum 
subspace of this system will evolves into a superposition 
of itself and the states in which the atom is unexcited 
and has released a photon into the CRW. With the half 
width of the band smaller than the driven frequency, the 
dynamic of the TLS is major governed by Hamiltonian 
in Eq.(fl"4]) with m — 0. In terms of the Green function, 
the probability for finding an initial excited TLS still in 
the excited state reads P e = 
we have defined 



dEe- %m C e (E) . Here 



C e (E) = (e<j> \(E-H R ) 



(16) 



which is the Fourier-Laplace transform of the amplitude 
C e (t) for the TLS in its excited state at arbitary time. 
Since the number of the cavities in the CRW is larger, 
the coupling between the TLS and the CRW is small. 
Therfore, Hamiltonian Hrq can be regarded as the un- 
perturbed part, while Hamiltonian Hr\ is treated as a 
perturbation part. Via Dyson's equation, we compute 
the Floquet-Green's functions up to the second nonvan- 
ishing order in the coupling strength g/y/N. Then the 
amplitude reads 



C e {E) = 



1 



E-ES 



, ~9n(E) 

E-El 



(17) 



where g n (E) is the Fourier-Laplace transform of 

g n (E) = / dtg n (t) e -^-W+™)r el Et- v t (lg) 
Jo 

9 2 Jlix) 



E 



N E-Et' 



The inverse Fourier-Laplace transform of g (E) yields the 
memory function 



9n(t) = g -J? l (x)Y, eU2 ' COSk 



(19) 



or reservoir response function [3J, |35j|, which depends on 
the quasiexcitation in the N modes of the CRW and char- 
acterizes the spectrum of the reservior. Comparing with 
the memory function $ (t) in the absence of modulation, 
an extra factor (x) has been involved. Obviously, by 
setting n = and A = 0, g n (t) reduces to $ (t). 

The inverse Fourier-Laplace transform of Eq. (|17[) 
yields the time evolution of the amplitude C e (t) 



C e = e l *° 



E't 



1-t I ,h[ \ -)g n (r)t ■ 



where detuning A = uj c — SI. The probability for finding 
the TLS in its excited state reads 



P e ~ exp (-Rt) . 



(21) 



Here the decay rate is the overlap of the modulation spec- 
trum /„ (w) and the reservoir coupling spectrum g n (w) 



R = 2tt 



dujf n (u>) g n (u>) , 



(22) 



where /„ (id) and g n (u>) is the Fourier transform of func- 
tions f n (r) and g n (r) . The function /„ (t) is defined 
as 



fn (r) = (l - -) e -^ A+n ^Q (t 



(23) 



where O (x) is the Heaviside unit step function, i.e., 
9 (x) = 1 for x > 0, and 9 (x) = for x < 0. Comparing 
with the form factor induced by the frequently measure- 
ment, an extra factor e~ mut has been introduced by the 
modulation. However one can find that when the mod- 
ulation is absent (n — 0), the modulation spectrum re- 
duces to the measurement-induced level-broadening func- 
tion in Ref.Q. 

The expression of functions /„ (t) in Eq. (|23p and g n (t) 
in Eq. p^l allows us to calculate the decay rate as 



R 



N 



Jn(x) 



k 



sin c' 



( A - 2£ cos k + nv) t 



(24) 



where sine = sin x jx. It shows that the decay rate R is 
determined by: 1) the parameters of the driving field, i.e. 
the driven intensity A and frequency v; 2) the detuning 
A between the cavity and the TLS; 3) the modulation 
time t; 4) the number N of cavities in the ID waveguide. 
But only frequency v, intensity A, and the detuning A 
can be adjusted experimentally. 

We first consider the long-time dynamics of the peri- 
odically driven TLS. As time t — > oo, the decay rate in 
Eq. (|2~41 is a sum of Dirac delta functions 



k 



2£ cos k + nv) 



(25) 



Eq. (125[) shows that depending on whether the matching 
condition 



A - 2£ cos k ■ 







(26) 
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is satisfied, the TLS is either (i) frozen to its initial ex- 
cited state or (ii) the TLS moves to the lower state and 
stays in it for ever. Case (ii) appears at the matching 
condition, but case (i) occurs when the transition energy 
of the TLS is out of resonance with the nth energy band 
of the CRW. When the number of cavities in the CRW 
is infinity, the states in the reservoir are continuum, one 
can replace the sums over k by integrals. Therefore, the 
properties of CRW are described by the reservoir spectral 
density 

p{w) = N- 1 ^ S(uj + 2^ cos k) (27) 

fe 

r o 2£<m 
= j oo 2£ = m 



And the Fourier transformation of the reservoir response 
function reads 

g n (uj)=g 2 jZ(A/v)p(u). (28) 

The behavior of the TLS is determined by whether the 
transition energy of the TLS is inside or outside the en- 
ergy band of the CRW. When Ek=o < — nv < £fc =7r , 
the TLS is in its ground state and the single quantum 
stays in the modes of the CRW. The decay rate reads 
R = 2-Kg n (A + nv). It is the extension of the golden rule 
rate to the case of a time-dependent coupling. When the 
TLS is out of resonance with the CRW, the TLS remains 
in its excited state. However, all the above matching 
conditions doesn't matter when the ratio of intensity A 
to frequency v is chosen such that J n {A/v) = 0. And 
at those points where the driving intensity A and fre- 
quency v satisfy J n {Ajv) — 0, the decay is completely 
suppressed. Actually, the occurrence of the complete sup- 
pression at J n {Ajv) = is caused by the decoupling be- 
tween the TLS and the CRW. Notice that the subscript 
n of the Bessel function is determined by the ratio of the 
energy difference between the CRW and the TLS to the 
frequency v, and the zeroes of the Bessel function J n (x) 
with different n appear in different argument x. There- 
fore, one can switch on or off the coupling between the 
TLS and the CRW by adjusting the detuning A for a 
given intensity A and v which satisfy J n {A/v) = 0. 

We now consider the dynamics of the periodically 
driven TLS with a finite time. The investigation starts 
with Eq.(f24|). Obviously, the decoupling between the 
TLS and the CRW still happens when the ratio of the 
driven intensity A to frequency v meets a root of the 
Bessel function J n (x). This decoupling preserves the 
population of the TLS. Consequently, QZE appears in 
this system. For J„ {A/v) ^ 0, one can still have QZE 
because the decay rate R grows as t increases when equa- 
tion (|26p is satisfied. However, When the transition en- 
ergy is out of resonance with the energy band of the 
CRW for any k, i.e. equation (j2"o) is not satisfied for 
any k, the decay rate is roughly a descending function 
of t. Consequently, the AZE occurs. Hence, for a given 





5 10 15 20 25 30 



FIG. 3. (Color online) The decay rate of the TLS as a function 
of time t with t in units of £ . The coupling strength g — 0.25 
and the number of cavities N = 41, A = l,A/v = 1 for 
blue dashed line, A = 3, A/v = 1 for red dot-dashed line, 
A = 3, A/v = 2.4 for green solid line. All parameters are in 
units of f. 

driven source, the Zeno-Anti-Zeno crossover can be ad- 
justed by the detuning between the TLS and the CRW. 
For a given detuning A, one can switch the occurrence of 
the QZE and AZE by controlling the driven intensity and 
frequency. In Fig|3J we plot the decay rate as a function 
of time t with the coupling strength g = 0.25£ and the 
number of the resonators N = 41. For the green solid 
line, the detuning A = 3, the ratio of the intensity to 
frequency A/v = 2.4. which satisfies the Jq(A/v) = 0, 
therefore the decay rate is vanished. For the red dot- 
dashed line, the transition energy of the TLS is outside 
the band of the CRW, hence, anti-Zeno effect appears. 
For the blue dashed line, transition energy inside the 
band, consequently, quantum Zeno effect occurs. There- 
fore, one can switch quantum Zeno effect to Anti-Zeno 
effect by tuning the transition energy of the TLS or the 
intensity and frequency of the driven field. 

As for the case with infinite cavities and finite time, 
the important factors are the the width and center of 
the spectrums, which determine the appearance of Zeno 
and anti-Zeno effect. We denote th width and center 
of the modulation spectrum as Af = t^ 1 and ojf = 
A + nv. The width and center of the reservoir cou- 
pling spectrum g n (ui) read A g = \[2^gJ n {A/v) and 
ujg = 0, respectively. When A/ 3> A g ,ujf, the effec- 
tive decay rate R ~ 27r/„ (A + nv), which grows with 
t, consequently, the quantum Zeno effect generally oc- 
curs. When Af <C \uf —u g \ and w/ is significantly de- 
tuned from the center oj g of g n (lo), the effective decay 
R ~ 27r<7„ (A + nv) <C 2irg n (oJ g ). As a result, the effec- 
tive decay R is a increasing function of the width Af, 
which leads to the acceleration of decay, i.e. the anti- 
Zeno effect. 



V. CONCLUSION 

In summary, we have considered a CRW with a period- 
driven TLS inside one of the resonators. The spontaneous 
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emission of the TLS is investigated in the one-quantum 
subspace via the Floquet representation of this system. 
Via computing the Floquet-Green's function up to the 
second nonvanishing order in the coupling strength, the 
amplitude for the TLS still in its excited state is obtained. 
It is found that the band structure of the CRW in the 
dispersion relation of the electromagnetic held gives rise 
to crossover of the quantum Zeno effect and anti-Zeno 
effect, which can be controlled by the driven intensity 
A and frequency v, the detuning A between the cavity 
and the TLS. One important effect induced by the diag- 
onal couplings between the driving held and the TLS is 
the strong dynamical suppression of decaying at suitable 
values of applied ac held. 
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NSFC No. 10775048, No. 10704023, and No. 11074071, 
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Project-sponsored by SRF for ROCS, SEM [2010]609- 
5, the Key Project of Chinese Ministry of Education 
(No. 210150), and Scientific Research Fund of Hu- 
nan Provincial Education Department No. 09B063, and 
No. 09C638. 

Appendix A: Floquet formulation 

Here, we give a brief outline of this method. Accord- 
ing to Floquet's theorem (4ll. l4ol |44]| . the time-dependent 
Schrodinger equation exists Floquet-state solutions of the 
form 

|*a (*)) = e- fe -* |$ a (*)) , (Al) 

where the Floquet state |$ Q (i)) = |$ Q (t + T)) is a func- 
tion with the same period as the driving field, e a is a 
real- valued energy function and is called the quasi-energy. 
The Floquet states obey the eigenvalue equation 

(h i^j \$ a (<)) = s a |* a (t)> . (A2) 

The operator on the left hand of the above equation gives 
the Floquet Hamiltonian 

H F = H-i— (A3) 
at 

which is time-independent. Floquet Hamiltonian oper- 
ates on an extended Hilbert space TZ = T-L ® T, which 
made up of the Hilbert space % of the system and the 



temporal space T of time-periodic functions. The tempo- 
ral part can be spanned by the orthonormal set of func- 
tions (t \m) = exp {imvt), where m = 0, ±1, ±2, ■ • • is the 
Fourier index, and 

1 f T 

{n\m} = — exp [i (m - n) vt] dt = S mn . (A4) 
1 Jo 

The composite Hilbert space TZ is called Floquet or 
Sambe space[43j. In Floquet space, the Floquet Hamil- 
tonian is linear and Hermitian, and the Floquet states 
provide a complete basis with the scalar product defined 
as 

«* Q (t) | 3> {£))) = i jT (4> a (*) ^ (*)) dt. (A5) 

The matrix elements of the time-evolution operator 
U a /3 (t, to) propagates the state \a) at time to to the state 
\(3) at time t > to according to the time-dependent Hamil- 
tonian H. In the Floquet representation, U a f}(t,ta) is 
related to the Floquet Hamiltonian 

U afi (t, t ) = (H ex P HHf (t - to)] |o0> e mwt 

n 

_ (A6) 
and is interpreted in Ref.p3|| as "the amplitude that 
a system initially in the Floquet state \a0) at time to 
evolves to the Floquet state \(3n) at time t according to 
the time-independent Floquet Hamiltonian Hp , summed 
over n with weighting factors exp (inujt)" . In experiment, 
the probability to go from the initial state |a) to the fi- 
nal state |/3) is the time- averaged transition probability 
between |a) and \f3) 

P a p = \^ n \ ex P H#f (t - t )} \a0)\ 2 . (A7) 

n 

The matrix elements U a p (t, to) is related to the Floquet- 
Green's functions via the Cauchy integral formula 
UaB (t, tp) = § e~ tEt G a pdE with the Floquet-Green's 
functions [4a] 

G^ = J]e tart Gg. (A8) 

n 

The right hand of equation (|A8|) is the Fourier expansion 
of G a p with the Fourier coefficients 

G [ 2\ = {H{E-H F )- l \aQ). (A9) 
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